Nemabiome metabarcoding shows a high prevalence of Haemonchus contortus and predominance of Camelostrongylus mentulatus in alpaca herds in the northern UK.

Gastrointestinal nematodes (GINs) are a common threat faced by pastoral livestock. Since their major introduction to the UK in the early 1990s, South American camelids have been cograzed with sheep, horses, and other livestock, allowing exposure to a range of GIN species. However, there have been no molecular-based studies to investigate the GIN populations present in these camelids. In the current study, we sampled nine alpaca herds from northern England and southern Scotland and used high-throughput metabarcoded sequencing to describe their GIN species composition. A total of 71 amplicon sequence variants (ASVs) were identified representing eight known GIN species. Haemonchus contortus was the most prevalent species found in almost all herds in significant proportions. The identification of H. contortus in other livestock species is unusual in the northern UK, implying that alpacas may be suitable hosts and potential reservoirs for infection in other hosts. In addition, the camelid-adapted GIN species Camelostrongylus mentulatus was identified predominantly in herds with higher faecal egg counts. These findings highlight the value of applying advanced molecular methods, such as nemabiome metabarcoding to describe the dynamics of gastrointestinal nematode infections in novel situations. The results provide a strong base for further studies involving cograzing animals to confirm the potential role of alpacas in transmitting GIN species between hosts.

Validation of deep amplicon sequencing of Dicrocoelium in small ruminants from Northern regions of Pakistan.

Dicrocoelium lancet flukes cause significant production loss in ruminant livestock. Although co-infection with multiple Dicrocoelium species within a host is common, techniques for studying the composition of these complex parasite communities are lacking. The pathogenicity, epidemiology, and therapeutic susceptibility of different helminth species vary, and little is known about the interactions that take place between co-infecting species and their hosts. Here, we describe the first applicationof metabarcoding deep amplicon sequencing method to studythe Dicrocoelium species in sheep and goats. First, rDNA ITS-2 sequences of four Dicrocoelium species (Dicrocoelium dendriticum, Dicrocoelium hospes, Dicrocoelium orientalis, and Dicrocoelium chinensis) were extracted from the NCBI public database. Phylogenetic analysis revealed separate clades of Dicrocoelium species; hence, molecular differentiation between each species is possible in co-infections. Second, 202 flukes belonging to seventeen host populations (morphologically verified as belonging to the Dicrocoelium genus) were evaluated to determine the deep amplicon sequencing read threshold of an individual fluke for each of the four species. The accuracy of the method in proportional quantification of samples collected from single hosts was further assessed. Overall, 198 (98.01%) flukes were confirmed as D. dendriticum and 1.98% produced no reads. The comparison of genetic distances between rDNA ITS-2 revealed 86% to 98% identity between the Dicrocoelium species. Phylogenetic analysis demonstrated a distinct clustering of species, apart from D. orientalis and D. chinensis, which sit very close to each other in a single large clade whereas D. hospes and D. dendriticum are separated into their own clade. In conclusion each sample was identified as D. dendriticum based on the proportion of MiSeq reads and validated the presence of this group of parasites in the Gilgit Baltistan and Khyber Pakhtunkhwa provinces of Pakistan. The metabarcoding deep amplicon sequencing technology and bioinformatics pathway have several potential applications, including species interactions during co-infections, identifying the host and geographical distribution of Dicrocoelium in livestock, drug therapy response evaluation and understanding of the emergence and spread of drug resistance.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds.

BACKGROUND: The genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome. RESULTS: We identified several thousand variants that had significantly steep clines ('SCV') across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle. CONCLUSIONS: To our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

A new record of the occurrence of Trichuris skrjabini Baskakov, 1924 in goats of Pakistan.

More than 23 Trichuroidea species have been identified in ruminants in different parts of the world. Most are pathogenic, causing trichurosis. Trichuris adults of most species within this family have a predilection for the ceca, where they may cause damage to the epithelial wall. In the present study, Trichuris spp. from large intestine of goats were analysed based on morphological and molecular characteristics. Fifty adult worms (male = 25 and female = 25) were selected for morphometric and molecular analysis. Male Trichuris were distinguished by their longer spicules, typical spicule sheaths, and small spicules that were always completely covered by the spicule sheath. The presence of an uneverted vulva in the vagina distinguished female worms. We have performed the molecular characterisation of adult warms to identify as Trichuris skrjabini. Genetic comparison of T. skrjabini rDNA ITS2 sequences with those from other Trichuris spp. was performed to assess within and between species variation and validate the use of ITS-2 rDNA as a robust species-specific marker for T. skrjabini identification. This work provides the first report of this parasite species from Pakistan and validated species-specific marker of T. skrjabini that reduces the production potential of goats in the country.

Histopathology and antibody responses describe the seasonal pattern of dicrocoeliosis in small ruminants in the Himalayan ranges of Pakistan.

In some parts of the world, Dicrocoelium spp. lancet flukes cause significant production loss in pastoral livestock, and accurate diagnosis of infection is important. The aims of the present study were to describe the histopathology and to investigate the transmission patterns of Dicrocoelium amongst ten sheep and goat farms in Khyber Pakhtunkhwa and Gilgit Baltistan, Pakistan. The liver histology and indirect enzyme-linked immunosorbent assay (ELISA) analyses followed standard procedures. The liver histopathology showed intensive tissue destruction and biliary hyperplasia associated with presence of adult flukes, severe inflammatory cell infiltration, congestion of blood vessels, damaged hepatocytes, and sinusoids in the infected areas. The time of onset of infection was investigated by ELISA detection of antibodies in sheep (n = 164) and goats (n = 152). Colostral transfer of Dicrocoelium antibodies from seropositive mothers was detected in sheep and goats up to 16 weeks of age. In both sheep and goats, the estimated time of infection differed between farms and years. Infection was seen in both sheep flocks and goat herds, with high variation between flocks and herds, and the highest infection rate in lambs. Dicrocoelium infection was most prevalent in sheep and goats in September (n = 84) and August (n = 63) respectively. This study concluded Dicrocoelium causes severe inflammation and necrosis of liver tissues in sheep and goats. Colostral transfer of antibodies can be detected up to about ten weeks of age. Higher infection rates are observed during August and September in sheep than in goats, putatively due to effects of different grazing and browsing behaviors on the ingestion of ants. The results will aid in the development of effective disease control strategies to ensure optimal growth and productivity of sheep and goats.

Spatial Distribution of Dicrocoelium in the Himalayan Ranges: Potential Impacts of Ecological Niches and Climatic Variables.

PURPOSE: Dicrocoeliosis can be an important cause of production loss in ruminants due to the cost of liver condemnation at slaughter. The aim of the present study was to determine the prevalence of Dicrocoelium infection and to predict the ecological niches and climatic variables that support dicrocoeliosis in the Himalayan ranges of Pakistan. METHODS AND RESULTS: Dicrocoelium was detected in 33 of 381 liver samples and 238 of 6060 blood samples taken from sheep and goat herds in the area. The prevalence of dicrocoeliosis was higher in sheep than in goats and highest in females aged more than 3 years. An environmental risk map was created to predict active zones of transmission and showed the highest probability values in central parts of the Chitral district in the northwest of Pakistan. Climatic variables of the mean monthly diurnal temperature range (Bio2), annual precipitation (Bio12), and normalised difference vegetation index (NDVI) were found to be significantly (p < 0.05) associated with the presence of Dicrocoelium infection. CONCLUSION: Together, the findings of this study demonstrate the most suitable ecological niches and climatic variables influencing the risk of dicrocoeliosis in the Himalayan ranges of Pakistan. The methods and results could be used as a reference to inform the control of dicrocoeliosis in the region.

Extensive testing of a multi-locus sequence typing scheme for Giardia duodenalis assemblage A confirms its good discriminatory power.

BACKGROUND: The flagellated parasite Giardia duodenalis is a major and global cause of diarrhoeal disease. Eight genetically very distinct groups, known as assemblages A to H, have been recognized in the G. duodenalis species complex, two of which (assemblages A and B) infect humans and other mammalian hosts. Informative typing schemes are essential to understand transmission pathways, characterize outbreaks and trace zoonotic transmission. In this study, we evaluated a published multi-locus sequence typing (MLST) scheme for G. duodenalis assemblage A, which is based on six polymorphic markers. METHODS: We genotyped 60 human-derived and 11 animal-derived G. duodenalis isolates collected in Europe and on other continents based on the published protocol. After retrieving previously published genotyping data and excluding isolates whose sequences showed allelic sequence heterozygosity, we analysed a dataset comprising 146 isolates. RESULTS: We identified novel variants at five of the six markers and identified 78 distinct MLST types in the overall dataset. Phylogenetic interpretation of typing data confirmed that sub-assemblage AII only comprises human-derived isolates, whereas sub-assemblage AI comprises all animal-derived isolates and a few human-derived isolates, suggesting limited zoonotic transmission. Within sub-assemblage AII, isolates from two outbreaks, which occurred in Sweden and Italy, respectively, had unique and distinct MLST types. Population genetic analysis showed a lack of clustering by geographical origin of the isolates. CONCLUSION: The MLST scheme evaluated provides sufficient discriminatory power for epidemiological studies of G. duodenalis assemblage A.

Faecal egg counts and nemabiome metabarcoding highlight the genomic complexity of equine cyathostomin communities and provide insight into their dynamics in a Scottish native pony herd.

Understanding the composition of gastrointestinal nematode communities may help to mitigate or exploit parasite adaptations within their host. We have used nemabiome deep amplicon sequencing of internal transcribed spacer-2 (ITS-2) ribosomal DNA to describe the temporal and host species composition of gastrointestinal nematode communities following sampling of six Scottish ponies across 57 months. In the absence of parasite control, each horse showed seasonal trends of increases and decreases in faecal egg counts, consistent with the epidemiology of equine strongylid parasites, however, the composition of parasites within individuals changed over time. Sixteen presumptive strongylid species were identified in each of the horses, 13 of which were distributed in a complex clade together with small numbers of amplicon sequences which could not be classified beyond the Cyathostominae subfamily level. Egg shedding of seven trichostrongylid species, which had previously been identified in co-grazed Soay sheep, was identified during the early spring. Faecal egg counts and the percentage of amplicon sequences assigned to each gastrointestinal nematode species were combined to describe their relative abundance across both host and time. Significant differences in species diversity between horses and between months were observed, being greatest from March to May and least from October to December. The magnitude of the individual horse effect varied between months and, conversely, the magnitude of the seasonal effect varied between individual horses. The most abundant gastrointestinal nematode in each of the horses was Cylicostephanus longibursatus (46.6% overall), while the abundance of the other strongylid species varied between horses and relative to each other. Patent C. longibursatus infections over the winter months might represent a genetic adaptation towards longer adult worm survival, or a lower rate of developmental arrest in the autumn. This study provides insight into highly complex phylogenetic relationships between closely related cyathostomin species; and describes the dynamics of egg shedding and pasture contamination of co-infecting equine gastrointestinal nematode communities. The results could be applied to determine how climatic and management factors affect the equilibrium between hosts and their parasites, and to inform the development of sustainable gastrointestinal nematode control strategies for different host species.

Genomic landscape of drug response reveals mediators of anthelmintic resistance.

Like other pathogens, parasitic helminths can rapidly evolve resistance to drug treatment. Understanding the genetic basis of anthelmintic drug resistance in parasitic nematodes is key to tracking its spread and improving the efficacy and sustainability of parasite control. Here, we use an in vivo genetic cross between drug-susceptible and multi-drug-resistant strains of Haemonchus contortus in a natural host-parasite system to simultaneously map resistance loci for the three major classes of anthelmintics. This approach identifies new alleles for resistance to benzimidazoles and levamisole and implicates the transcription factor cky-1 in ivermectin resistance. This gene is within a locus under selection in ivermectin-resistant populations worldwide; expression analyses and functional validation using knockdown experiments support that cky-1 is associated with ivermectin survival. Our work demonstrates the feasibility of high-resolution forward genetics in a parasitic nematode and identifies variants for the development of molecular diagnostics to combat drug resistance in the field.

Genetic characterisation of the Theileria annulata cytochrome b locus and its impact on buparvaquone resistance in bovine.

Control of tropical theileriosis, caused by the apicomplexan Theileria annulata, depends on the use of a single drug, buparvaquone, the efficacy of which is compromised by the emergence of resistance. The present study was undertaken to improve understanding of the role of mutations conferring buparvaquone resistance in T. annulata, and the effects of selection pressures on their emergence and spread. First, we investigated genetic characteristics of the cytochrome b locus associated with buparvaquone resistance in 10 susceptible and 7 resistant T. annulata isolates. The 129G (GGC) mutation was found in the Q01 binding pocket and 253S (TCT) and 262S (TCA) mutations were identified within the Q02 binding pocket. Next, we examined field isolates and identified cytochrome b mutations 129G (GGC), 253S (TCT) and 262S (TCA) in 21/75 buffalo-derived and 19/119 cattle-derived T. annulata isolates, providing evidence of positive selection pressure. Both hard and soft selective sweeps were identified, with striking differences between isolates. For example, 19 buffalo-derived and 7 cattle-derived isolates contained 129G (GGC) and 253S (TCT) resistance haplotypes at a high frequency, implying the emergence of resistance by a single mutation. Two buffalo-derived and 12 cattle-derived isolates contained equally high frequencies of 129G (GGC), 253S (TCT), 129G (GGC)/253S (TCT) and 262S (TCA) resistance haplotypes, implying the emergence of resistance by pre-existing or recurrent mutations. Phylogenetic analysis further revealed that 9 and 21 unique haplotypes in buffalo and cattle-derived isolates were present in a single lineage, suggesting a single origin. We propose that animal migration between farms is an important factor in the spread of buparvaquone resistance in endemic regions of Pakistan. The overall outcomes will be useful in understanding how drug resistance emerges and spreads, and this information will help design strategies to optimise the use and lifespan of the single most drug use to control tropical theileriosis.

Genomic signatures of selection associated with benzimidazole drug treatments in Haemonchus contortus field populations.

Genome-wide methods offer a powerful approach to detect signatures of drug selection. However, limited availability of suitable reference genomes and the difficulty of obtaining field populations with well-defined, distinct drug treatment histories mean there is little information on the signatures of selection in parasitic nematodes and on how best to detect them. This study addresses these knowledge gaps by using field populations of Haemonchus contortus with well-defined benzimidazole treatment histories, leveraging a recently completed chromosomal-scale reference genome assembly. We generated a panel of 49,393 genomic markers to genotype 20 individual adult worms from each of four H. contortus populations: two from closed sheep flocks with an approximate 20 year history of frequent benzimidazole treatment, and two populations with a history of little or no treatment. Sampling occurred in the same geographical region to limit genetic differentiation and maximise the detection sensitivity. A clear signature of selection was detected on chromosome I, centred on the isotype-1 ß-tubulin gene. Two additional, but weaker, signatures of selection were detected; one near the middle of chromosome I spanning 3.75 Mbp and 259 annotated genes, and one on chromosome II spanning a region of 3.3 Mbp and 206 annotated genes, including the isotype-2 ß-tubulin locus. We also assessed how sensitivity was impacted by sequencing depth, worm number, and pooled versus individual worm sequence data. This study provides the first known direct genome-wide evidence for any parasitic nematode, that the isotype-1 ß-tubulin gene is quantitatively the single most important benzimidazole resistance locus. It also identified two additional genomic regions that likely contain benzimidazole resistance loci of secondary importance. This study provides an experimental framework to maximise the power of genome-wide approaches to detect signatures of selection driven by anthelmintic drug treatments in field populations of parasitic nematodes.

Proteomic Comparison of Ivermectin Sensitive and Resistant Staphylococcus aureus Clinical Isolates Reveals Key Efflux Pumps as Possible Resistance Determinants.

Ivermectin (IVM) is a versatile drug used against many microorganisms. Staphylococcus aureus is one of the most devastating microorganisms. IVM sensitive and resistant S. aureus strains were recently reported. However, the underlying molecular mechanisms of resistance are unknown. Clinical isolates of S. aureus were used for determination of the sensitivities against IVM by growth curve analysis and time-kill kinetics. Then, proteomic, and biochemical approaches were applied to investigate the possible mechanisms of resistance. Proteomic results showed a total of 1849 proteins in the dataset for both strains, 425 unique proteins in strain O9 (IVM sensitive), and 354 unique proteins in strain O20 (IVM resistant). Eight proteins with transport functions were differentially expressed in the IVM resistant strain. Among them, three efflux pumps (mepA, emrB, and swrC) were confirmed by qPCR. The IVM resistant S. aureus may overexpress these proteins as a key resistance determinant. Further experiments are required to confirm the exact mechanistic relationship. Nevertheless, the possibility of blocking these transporters to reverse or delay the onset of resistance and reduce selection pressure is potentially appealing.

Transcriptomic analyses implicate neuronal plasticity and chloride homeostasis in ivermectin resistance and response to treatment in a parasitic nematode.

The antiparasitic drug ivermectin plays an essential role in human and animal health globally. However, ivermectin resistance is widespread in veterinary helminths and there are growing concerns of sub-optimal responses to treatment in related helminths of humans. Despite decades of research, the genetic mechanisms underlying ivermectin resistance are poorly understood in parasitic helminths. This reflects significant uncertainty regarding the mode of action of ivermectin in parasitic helminths, and the genetic complexity of these organisms; parasitic helminths have large, rapidly evolving genomes and differences in evolutionary history and genetic background can confound comparisons between resistant and susceptible populations. We undertook a controlled genetic cross of a multi-drug resistant and a susceptible reference isolate of Haemonchus contortus, an economically important gastrointestinal nematode of sheep, and ivermectin-selected the F2 population for comparison with an untreated F2 control. RNA-seq analyses of male and female adults of all populations identified high transcriptomic differentiation between parental isolates, which was significantly reduced in the F2, allowing differences associated specifically with ivermectin resistance to be identified. In all resistant populations, there was constitutive upregulation of a single gene, HCON_00155390:cky-1, a putative pharyngeal-expressed transcription factor, in a narrow locus on chromosome V previously shown to be under ivermectin selection. In addition, we detected sex-specific differences in gene expression between resistant and susceptible populations, including constitutive upregulation of a P-glycoprotein, HCON_00162780:pgp-11, in resistant males only. After ivermectin selection, we identified differential expression of genes with roles in neuronal function and chloride homeostasis, which is consistent with an adaptive response to ivermectin-induced hyperpolarisation of neuromuscular cells. Overall, we show the utility of a genetic cross to identify differences in gene expression that are specific to ivermectin selection and provide a framework to better understand ivermectin resistance and response to treatment in parasitic helminths.

Current methods for the detection of Plasmodium parasite species infecting humans.

Malaria is the world's fatal parasitic disease. The ability to quickly and accurately identify malaria infection in challenging environments is crucial to allow efficient administration of the best treatment regime for human patients. If those techniques are accessible and efficient, global detection of Plasmodium species will become more sensitive, allowing faster and more precise action to be taken for disease control strategies. Recent advances in technology have enhanced our ability to diagnose different species of Plasmodium parasites with greater sensitivity and specificity. This literature review provides a summary and discussion of the current methods for the diagnosis and identification of Plasmodium spp. in human blood samples. So far not a single method is precise, but advanced technologies give consistent identification of a Plasmodium infection in endemic regions. By using the power of the recent methods, we can provide a broader understanding of the multiplicity of infection and or transmission dynamics of Plasmodium spp. This will result in improved disease control strategies, better-informed policy, and effective treatment for malaria-positive patients.

Molecular Characterization of Toxoplasma gondii in Cats and Its Zoonotic Potential for Public Health Significance.

Toxoplasmosis is a globally distributed disease of warm-blooded animals. It is caused by the opportunistic parasite Toxoplasma gondii (T. gondii). One-third of the global human population is believed to be infected with T. gondii. Cats serve as final host of T. gondii and are the main source of contamination of soil and water. This study aimed to detect genotypes of T. gondii in cats. Fecal samples (n = 400) were collected from districts of South Punjab (Khanewal and Sahiwal), and were processed by polymerase chain reaction (PCR) followed by sequencing and phylogenetic analysis. The obtained oligonucleotide sequences (T. gondii) were submitted to the GenBank database, and the evolutionary tree was constructed using MEGA-X software. Seven fecal samples (3.5%) from cats were positive. Five out of thirteen fecal samples (38.46%) found to be positive for T. gondii with microscopy were confirmed by PCR. After phylogenetic analysis with 3 clonal types and atypical strains, isolates of T. gondii in current study were more closely linked to a typical strain (AF249696). Besides genotyping from cats, seroprevalence from humans and ruminants is still considered to be the best and easiest way to identify the Toxoplasma. Blood samples were collected from sheep and goats (n = 2000 each), and human blood samples (n = 400) were collected from the same vicinity. Seroprevalence was determined using a commercial enzyme-linked immunosorbent assay (ELISA) kit. In Khanewal, the blood samples of 292 goats (29.2%) and 265 sheep (26.5%), and 6 fecal samples from cats (3%) were positive. Out of 200 human blood samples, 52 were positive, with a seroprevalence of 26%. In the Sahiwal district, the blood samples from 49 humans, 235 sheep and 348 goats were positive, with seroprevalence of 24.5%, 23.5% and 34.8%, respectively. The present study revealed the current circulating genotype of T. gondii from cats in the districts Khanewal and Sahiwal and the seroprevalence of the organism in small ruminants and humans living in the same vicinity. Further genotype analyses of the organism from ruminants and humans are needed.

Contamination of Soil, Water, Fresh Produce, and Bivalve Mollusks with Toxoplasma gondii Oocysts: A Systematic Review.

Toxoplasma gondii is a major foodborne pathogen capable of infecting all warm-blooded animals, including humans. Although oocyst-associated toxoplasmosis outbreaks have been documented, the relevance of the environmental transmission route remains poorly investigated. Thus, we carried out an extensive systematic review on T. gondii oocyst contamination of soil, water, fresh produce, and mollusk bivalves, following the PRISMA guidelines. Studies published up to the end of 2020 were searched for in public databases and screened. The reference sections of the selected articles were examined to identify additional studies. A total of 102 out of 3201 articles were selected: 34 articles focused on soil, 40 focused on water, 23 focused on fresh produce (vegetables/fruits), and 21 focused on bivalve mollusks. Toxoplasma gondii oocysts were found in all matrices worldwide, with detection rates ranging from 0.09% (1/1109) to 100% (8/8) using bioassay or PCR-based detection methods. There was a high heterogeneity (I2 = 98.9%), which was influenced by both the sampling strategy (e.g., sampling site and sample type, sample composition, sample origin, season, number of samples, cat presence) and methodology (recovery and detection methods). Harmonized approaches are needed for the detection of T. gondii in different environmental matrices in order to obtain robust and comparable results.

A novel metabarcoded deep amplicon sequencing tool for disease surveillance and determining the species composition of Trypanosoma in cattle and other farm animals.

The World Health Organization (WHO) and the Food and Agriculture Organization (FAO) have developed strategies to control trypanosomiasis in humans and livestock in endemic areas. These require a better understanding of the distribution of different Trypanosoma species and improved predictions of where they might appear in the future, based on accurate diagnosis and robust surveillance systems. Here, we describe a metabarcoding deep amplicon sequencing method to identify and determine the Trypanosoma species in co-infecting communities. First, four morphological verified Trypanosoma species (T. brucei, T. congolense, T. vivax and T. theileri) were used to prepare test DNA pools derived from different numbers of parasites to evaluate the method's detection threshold for each of the four species and to assess the accuracy of their proportional quantification. Having demonstrated the accurate determination of species composition in Trypanosoma communities, the method was applied to determine its detection threshold using blood samples collected from cattle with confirmed Trypanosoma infections based on a PCR assay. Each sample showed a different Trypanosoma species composition based on the proportion of MiSeq reads. Finally, we applied the assay to field samples to develop new insight into the species composition of Trypanosoma communities in cattle, camels, buffalo, horses, sheep, and goat in endemically infected regions of Pakistan. We confirmed that Trypanosoma evansi is the major species in Pakistan and for the first time showed the presence of Trypanosoma theileri. The metabarcoding deep amplicon sequencing method and bioinformatics pathway have several potential applications in animal and human research, including evaluation of drug treatment responses, understanding of the emergence and spread of drug resistance, and description of species interactions during co-infections and determination of host and geographic distribution of trypanosomiasis in humans and livestock.

A loop-mediated isothermal amplification (LAMP) assay to identify isotype 1 ß-tubulin locus SNPs in synthetic double-stranded Haemonchus contortus DNA.

Development of sustainable gastrointestinal nematode (GIN) control strategies depends on the ability to identify the frequencies of drug-susceptible and resistant genotypes in GIN populations arising from management practices undertaken on individual farms. Resistance to BZ drugs in GINs has been shown to be conferred by the presence of defined SNPs in the isotype 1 ß-tubulin locus. Loop-mediated isothermal amplification (LAMP) assays are amenable to use on a range of DNA templates and are potentially adaptable to use in practical, cost-effective, pen-side diagnostic platforms that are needed to detect anthelmintic resistance in the field. In this study, we designed primers and examined LAMP assays to detect each of the three major isotype 1 ß-tubulin SNPs conferring genetic susceptibility to BZ drugs. We used artificial pools of synthetic DNA, containing different proportions of susceptible and resistant SNPs to determine reproducibility of the assays. We demonstrated the detection of each of the isotype 1 ß-tubulin SNPs conferring susceptibility to BZ drugs using the optimal LAMP assay. Isotype 1 ß-tubulin SNP typing was effective in detecting BZ susceptibility, but the accuracy was reduced in samples with less than 60 % susceptible DNA. Our results show the potential for LAMP SNP typing to detect genetic susceptibility or resistance to anthelmintic drugs in livestock GINs, and some of the limitations in our approach that will need to be overcome in order to evaluate this assay using field samples. Supplementary Information: The online version contains supplementary material available at 10.1007/s12639-021-01414-w.

High-throughput sequencing of Fasciola spp. shows co-infection and intermediate forms in Balochistan, but only Fasciola gigantica in the Punjab province of Pakistan.

Fasciola gigantica and Fasciola hepatica are digenetic trematodes causing fasciolosis in ruminants. The host and geographical distribution of both Fasciola species are influenced by environmental and climatic conditions favouring survival and development of free-living stages and intermediate hosts, and livestock management practices. The aim of the present study was to describe the host distribution of the two Fasciola species in buffalo, cattle, goats, and sheep in the Balochistan and Punjab provinces of Pakistan. 359 flukes were collected from a total of 32 livers from the four livestock species. Deep amplicon sequencing of the internal transcribed spacer region 2 of ribosomal DNA (rDNA ITS-2) and mitochondrial nicotinamide adenine dinucleotide dehydrogenase 1 (mtDNA ND-1) loci confirmed co-infection of F. hepatica and F. gigantica in Balochistan and single species F. gigantica infection in Punjab. In Balochistan, co-infections and hybrids of both Fasciola species were identified in cattle, with more F. hepatica detected than F. gigantica. However, F. hepatica was the only species identified in goats, and F. gigantica was the only species identified in buffalo. In Punjab, all flukes were confirmed as F. gigantica in each of the four livestock species. Overall, the results indicate differences in the host and geographical distribution of F. gigantica and F. hepatica, and provide useful knowledge for the development of control strategies for livestock and humans.

Molecular characterization and phylogenetic analyses of Fasciola gigantica of buffaloes and goats in Punjab, Pakistan.

Fasciola gigantica is considered to be a major pathogen causing fasciolosis in the Indian subcontinent, resulting in production losses of millions of dollars in the livestock industry. Understading the dispersal origin and the patterns of spread of F. gigantica is important. A total of 53 Fasciola flukes collected from buffaloes and goats in Punjab, Pakistan between 2017 and 2018 were identified as F. gigantica based on the multiplex PCR for the phosphoenolpyruvate carboxykinase (pepck) and the PCR-restriction fragment length polymorphism (RFLP) for DNA polymerase delta (pold). A significant genetic difference between F. gigantica from buffaloes and goats was indicated by the genetic analyses of mitochondrial markers, NADH dehydrogenase subunit 1 (nad1) and cytochrome C oxidase subunit 1 (cox1). Phylogenetic analysis of the seventeen nad1 haplotypes of F. gigantica from Pakistan with those in neighbouring countries of the Indian subcontinent revealed that all the haplotypes identified in Pakistan were clustered in haplogroup A. fasciola gigantica with the eight haplotypes might be expanded in Pakistan from Indian origin, along with the migration of the domestic animals, since they were related to Indian haplotypes. In contrast, the remaining nine haplotypes were not shared with any neighbouring countries, suggesting independent origin, probably from neighbouring Middle East countries. However, cautious interpretation is required due to the very limited samples size of this study. Our study provides a proof of concept for a method that could be used to investigate the epidemiology of F. gigantica.